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ABSTRACT 

The faint end of the differential galaxy number counts, n(m), in the Hub- 
ble Deep Field (HDF) North has been determined for the F450W, F606W, and 
F814W filters by means of surface-brightness fluctuation (SBF) measurements. 
This technique allows us to explore n(m) beyond the limiting magnitude of the 
HDF, providing new, stronger constraints on the faint end of n(m). This has 
allowed us to test the validity of previous number count studies and to produce 
a new determination of the faint end of n(m) for magnitudes fainter than 28.8 
in the AB system and to extend this estimate down to 31. This value represents 
an extension of more than two magnitudes beyond the limits of previous photo- 
metric studies. The obtained n(m) slopes are 7 = 0.27, 0.21, and 0.26 in B i5Q , 
Vf306, and Isu, respectively. 

Subject headings: galaxies: evolution - galaxies: formation - cosmology: observa- 
tions 



1. Introduction 

The galaxy luminosity function ($) remains at the core of both galaxy evolution and 
cosmology. By integrating $ over space and time, various observable distributions can be 
obtained. In particular, differential number counts of galaxies as a function of apparent 
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magnitude, n(m), is obtained by integrating <£> over all redshifts and morphological types. 
The reliability of the predicted n{m) values depends directly on the cosmological model and 
on how well the evolution of galaxies is modelled from their formation to the present. The 
study of number counts can therefore be used to test world models or to search for evolution 
during the look-back time. 

The function n(m) provides one of the most fundamental observables and has been 
studied by several authors (for reviews see Ellis 1997; Koo & Kron 1992; Sandage 1988). 
Recently, a major effort has been made to reach very faint magnitudes using the Hubble 
Space Telescope (HST). Although the HST cannot compete with ground-based observa- 
tions in terms of collecting area, it does provide an unprecedent view of the optical sky at 
small angular scales and faint flux levels. A number of authors have studied faint galaxy 
counts based on the Hubble Deep Field (HDF) (Williams et al. 1996) and using different 
photometry packages (Williams et al. 1996; Metcalfe et al. 1996; Lanzzetta, Yahil, & 
Fernandez-Soto 1996; Pozzetti et al. 1998; Metcalfe et al. 2001). In this paper, we use 
the notation -B450, Vqqq, and ^814 to denote magnitudes in the HST passbands on the AB 
system (Oke 1974). Ferguson (1998) has compared the HDF catalogues available in the 
literature. They generally show good agreement, although the effects of different isophotal 
thresholds and different splitting algorithms are apparent. However, Ferguson, Dickinson, & 
Williams (2000) showed that, at J 8 i4 = 26, the different galaxy counts agree to within 25% 
in all the catalogues, whereas at I 8 u = 28 there is a factor of 1.7 difference among them. 
This emphasizes the fact that galaxy counting is not a precise science. 

In this paper, the HDF will be used to test the validity of previous number count studies 
and to produce a new determination of the faint end of n{m) for magnitudes fainter than 
28.8. Surface-brightness fluctuation (SBF) measurements are used. This allows us to explore 
n(m) beyond the limiting magnitude of the HDF and to overcome most of the limitations 
arising from incompleteness, providing new, stronger constraints on the faint end of n{m). 

2. The data 

This work is based on data from the HDF North. Observations were made in 1995 
December 18-30, and both raw and reduced data have been put into public domain as a 
community service (Williams et al. 1996). Version 2 F450W, F606W and F814W images, 
released into public domain on 1996 February 29, have been used here. 

The final Version 2 images were combined using the DRIZZLE algorithm. Drizzling 
causes the noise in one pixel to be correlated with the noise in the adjacent one (Williams 
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ct al. 1996). The SBF technique is based on the spectral analysis of an image. Because 
the power spectrum of the image is modified by the drizzling process, drizzled images were 
not considered in this study. "Weighted and cosmic-ray cleaned" images have been used 
instead. These are flat-fielded, cosmic-ray- rejected and sky-subtracted stacked images at 
each dither position. Only dark exposures have been selected. Each of the the three wide- 
field and the planetary camera WFPC2 chips were analyzed separately. The labels of all 
images considered, the corresponding filter and the total exposure times are listed in table 
1. 



3. Surface-Brightness fluctuations in the HDF-N 

The SBF concept was introduced by Tonry & Schneider (1988), who noted that, in 
the surface photometry of a galaxy far enough away to remain unresolved, a pixel-to-pixel 
fluctuation is observed because of the Poisson statistics of the spatial distribution of stars, 
globular clusters, background galaxies, etc. This technique was introduced with the aim of 
measuring distances. Comparing SBFs produced by the stellar population of a galaxy with 
those of nearby galaxies for which externally calibrated distances are available, accurate 
estimates of distances can be obtained up to ~40 Mpc (Tonry et al. 2000, 2001). Other 
authors have used SBF studies to determine the age and metallicity of unresolved stellar 
populations of nearby galaxies (Liu, Chariot, & Graham 2000; Blakeslee, Vazdekis, & Ajhar 
2001; Hidalgo, Marin- Franch, & Aparicio 2003); however, the SBF signal can provide in- 
formation about other kinds of undetected and unresolved objects in an image. This is the 
case, for example, of globular cluster populations (Blakeslee & Tonry 1995; Blakeslee 1999; 
Marin- Franch & Aparicio 2002, 2003). In this paper, SBFs have been used to characterize 
the faint end of n(m) using undetected galaxies in the HDF-N images. 

Next in this section, the theoretical background of SBFs and the HDF-N signal mea- 
surement are described in detail. 



3.1. Theory 

The SBF technique involves spectral analysis of the pixel-to-pixel fluctuation signal. 
This provides the total point spread function PSF-convolved variance (Po) produced by all 
point objects whose spatial flux distribution is convolved with the PSF, and the total non- 
PSF-convolved variance (-Pi). A convolution in the real space transforms into a product in 
Fourier space. For this reason, the power spectrum of an image, P(k), has the form: 
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P(k) = P 1 + P E(k), (1) 

where E(k) is the power spectrum of the PSF convolved with some window function. The 
E(k) can be approximated by the power spectrum of the PSF alone, Ppsf(^), with a negligible 
error (Jensen et al. 1998). So eq. 1 transforms into 



P(k) = P 1 + P P PSF (k). (2) 

Once the power spectrum of an image has been computed, the variances Pq and Pi can 
be obtained fitting eq. 2 to data. 

In HDF images, Pq is mainly produced by faint galaxies, and, as we will see later, by 
cosmic rays. So P must be equal to the sum of the variances produced by faint galaxies 
(cr| G ) and by cosmic rays (cr^ r ): 

Po = cr 2 BG + a 2 cr . (3) 

On the other hand, Pi is the sum of read-out noise (of ), photon shot noise (<7p h ), and dark 
current (cr^ c ) variances: 



Pi = a 2 Q + a 2 ph + a 2 dc . (4) 

The pixel-to-pixel variance produced by a class of objects can be evaluated as the second 
moment of the differential number counts of that object population (Tonry & Schneider 
1988). In the present study, the target population is composed by faint undetected galaxies. 
The brightest individuals, which are detected, must be masked out before the SBF analysis. 
If all sources brighter than a limiting flux (/u m ) are masked, then the variance produced by 
the remaining non-masked population is: 

= / <f)f 2 df. (5) 
Jo 

We can put this equation in terms of magnitudes via the relationship 

m= -2.51og(/) +m*, (6) 

/ being the flux (DN s _1 pix^ 1 ), and m* the magnitude of an object yielding one unit of 
flux per unit time; that is, the photometric zero point listed in Table 2. 

If n{m) is known, then the theoretical variance produced by faint galaxies can be esti- 
mated. Assuming the following pure power-law form for n{m): 

n(m) = A W m , (7) 
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where A is a normalizing constant and 7 is the slope of the magnitude distribution, the 
variance produced by the non-masked galaxy population is then 



a 



A 



BG 



In 10(0.8 -7) 



^Q7mc^g0.8(m*-m c )^ ^ 



where m c is the limiting magnitude; i.e. the magnitude corresponding to the limiting flux, 
/iim- We will refer to this variance as the n(m)-estimated <r| G from here on. As we have 
shown, in order to compute it, an initial n(m), obtained from a given photometric catalogue, 
must be assumed. 



3.1.1. Photometric catalogues 

A number of authors have studied n{m) in the HDF-N. The main discrepancies be- 
tween previous number count studies occur between Williams et al. (1996) and Metcalfe 
et al. (2001). These coincide for "bright" magnitudes, but large differences arise for fainter 
galaxies, the Metcalfe et al. (2001) number counts being larger than those of Williams et 
al. (1996). 

As pointed out by Ferguson (1998), there are two reasons for this discrepancy. First, 
as we go fainter the Metcalfe et al. (2001) magnitudes become systematically brighter than 
those of Williams et al. (1996). Metcalfe et al. (2001) claimed that the effect of this on the 
counts is actually quite small, generally < 10%. Second, Metcalfe et al. (2001) find objects 
that Williams et al. (1996) do not detect at all. This appears to account for the majority 
of the differences between the data sets. Metcalfe et al. (2001) argued that virtually all of 
these objects are merged in the reductions of Williams et al. (1996) but not in their data. 

In Figure 1, differential number counts results from Williams et al. (1996) (filled circles) 
and Metcalfe et al. (2001) (open circles) are plotted for the F450W, F606W, and F814W 
filters, respectively. These data have been obtained from tables 9 and 10 in Williams et 
al. (1996), and from tables 8, 10 and 12 in Metcalfe et al. (2001), all expressed in the 
total magnitude scale. Following eq. 7, n(m) has been fitted to the data of Williams et al. 
(1996) and Metcalfe et al. (2001), and the fitted functions (solid lines) are also plotted. As 
Williams et al. (1996) found a change in the n(m) slope at a magnitude of around 26, n(m) 
has been fitted to their data in the two magnitude intervals [23, 26] and [26, 29]. 

The fitted n{m) functions have been used in eq. 8 to compute the n(m)-estimated <r| G . 
Results are listed in table 3. These n(m)-estimated <x| G values will be later on compared with 
those directly derived from the SBF measurements (which we will call the SBF-measured 
<t| g ). This comparison will allow us to evaluate the validity of the differential number counts 
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of both Williams et al. (1996) and Metcalfe et al. (2001) and, as a result, a final n(m) will 
be proposed. 

Before describing the details of the SBF measurements in HDF images, note that the 
SBF technique is valid only if faint galaxies have a stellar appearance. In Ferguson (1998) 
the radius-magnitude relation for galaxies in the Williams et al. (1996) HDF catalogue 
has been analyzed. For magnitudes fainter than Vqog = 28.8, all the galaxies have a radius 
smaller that 0.16", close to the FWHM of the WF chips. So very faint galaxies, in the 
magnitude range where SBF will be measured, can be assumed to have a stellar appearance. 

3.2. Procedure to obtain SBF 

Here, the practical procedure for obtaining the SBF signal in HDF-N images is described 
in detail. First of all, it should be noted that cosmic rays are difficult to discriminate from 
stars in HST frames. In order to estimate and eliminate the cosmic-ray contribution to the 
SBF signal, Pq, a procedure based on the random nature of cosmic-ray events has been used 
for each filter and chip. The SBF signal has be measured not only on the final combined 
images, but on all the individual images listed in Table 1 as well. 

Before computing the power spectrum of an image, objects brighter than m c must be 
masked out. In this study, the window functions have been created using the Williams 
et al. (1996) photometric catalogue, which is the only one available to us. As isophotal 
magnitudes were considered while creating the window functions, in order to convert them to 
the total magnitude scale an isophotal-to-total magnitude correction of 0.2 mag (Williams et 
al. 1996) was applied. The SBF analysis have been be performed considering two different 
values of m c : 27.8 and 28.8. All objects brighter than m c have been masked out using a 
window function whose pixel values are zero in a circle centered on the location of the bright 
objects and unity in the rest of the image. The window function has been created using 
a patch radius large enough to completely mask bright galaxies, including their external 
haloes and therefore merged galaxies where these exist. The procedure creating the mask 
has been the following: first, the brightest galaxies have been masked one by one manually. 
In order to avoid residual light beyond the masked regions, very generous patch sizes have 
been adopted. The shape of the used patchs for these very bright galaxies depends on the 
shape of the particular masked galaxy. Once bright galaxies have been masked, the rest of 
galaxies brighter than the chosen m c have been masked using circular patches. The radius 
of these patches has been chosen to be the same for all galaxies, and its size is again very 
generous: the adopted patch radius is 15 pixels. 
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In order to test if all the residual light beyond the masked regions has been eliminated, 
the SBF analysis has been repeated for one image varying the radius of the patches. We 
have considered patch radii of 20, 25 and 30 pixels. Note that if a radius larger than 30 
pixels would be used, the image would be completely masked due to the superposition of 
adjacent patches. The considered image has been the F450W average image of WF2 with 
m c = 27.8. The SBF results are listed in table 4. It can be seen that the SBF results are 
independent of the patch radius. As a conclusion, it can be seen that the adopted patch 
radius avoid residual light beyond the masked regions, as required. 

If the differences between the Williams et al. (1996) and Metcalfe et al. (2001) data 
sets relies on the objects that are merged in the former and not in the latter, as argued by 
Metcalfe et al. (2001), then the created window function is virtually the same as one would 
obtain using the Metcalfe et al. (2001) photometric catalogue, and also masks objects that 
are in this catalogue and not in that of Williams et al. (1996). 

Next, multiplying the window function by the image, the residual masked image is 
obtained. It is this image which is used to compute the power spectrum. As has been said, 
two sets of masked images have been computed, one for each m c value. 

The power spectrum of the masked image is two-dimensional. It is radially averaged 
in order to obtain the one-dimensional power spectrum. Fitting the power spectrum of the 
images with eq. 2, the quantities P an d -Pi can be obtained for each image. As we could 
not construct a PSF from HDF images, since there were not enough stars in any of the four 
chips, the used PSFs were the high-S/N PSFs extracted by P. B. Stetson from a large set of 
uncrowded and unsaturated WFPC2 images. 

Spatial variations of the PSF along the CCD mosaic could introduce a significant uncer- 
tainty in P . To limit this effect, we have used a PSF template for each one of the four chips 
of the WFPC2. Minor PSF variations inside each particular chip have not been considered. 
For each chip, both PSF and the computed power spectrum represent mean values across 
the complete field of the chip, so any spatial variation of the PSF across the chip would 
affect the SBF fitting procedure, introducing an uncertainty in the Pq measurement. This 
uncertainty is small, and in any case, it is included in the obtained Pq uncertainty. 

In Fig. 2 we show an example of the SBF fitting procedure in an HDF image. The 
observed discrepancy at low wave numbers between the obtained power spectrum and the 
fit is due to large scale fluctuations in the background brightness of the images. This wave 
number region, where the discrepancy occurs, is not taken into account when fitting eq. 2. 
The power-spectrum fitting procedure is the following: eq. 2 is used to fit the power spectrum 
for wave numbers in the range [ko, k max ], /c max being the highest wave number of the computed 
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power spectrum, and k a number which we vary from to k m£LX . As a result, two functions, 
Po(k ) and Pi(k ), are obtained. The function -Po(^o) is also shown in Fig. 2 (small boxes). 
It can be seen that this function exhibits a "plateau" region. The final adopted result and 
its uncertainty for P are obtained computing the average and standard deviation of Po(^o) 
in the plateau interval. For the Pi measurement, the procedure is exactly the same as for 
Po- 

Pq and Pi results corresponding to m c = 27.8 and m c = 28.8 are listed in Tables 5 and 
6, respectively. Results are listed for all images, corresponding to all chips and all filters. 

4. Results 

In this section, results from Tables 5 and 6 will first be carefully analyzed, and will be 
used to test the probable effect of flat-fielding errors on the measured Po. Then, the obtained 
Po values in different images will be used to estimate the contribution of cosmic rays and 
to deduce the desired SBF- measured cr BG . Finally, in order to check the SBF results, two 
consistency tests will be performed. 

4.1. The effect of flat-fielding errors on P 

The possibility of flat-fielding errors affecting the measured Po must considered. If 
flat-fielding would contribute to Po, then its effect should be larger in images with high 
sky background. In this context, considering same filter and chip images, if flat-fielding is 
contributing to Po, a relation between Po and the sky level should appear. 

To test if the flat-fielding errors have an influence on the measured Po, filter F606W has 
been considered because it provides a larger number of images covering a wider range of sky 
levels. Sky levels are provided in the weighted, cosmic-ray cleaned image headers. As an 
example, fig. 3 shows the measured Po as a function of these sky levels for the WF2 F606W 
images with m c = 27.8. It can be seen that no relation is apparent, so flat-fielding errors are 
insignificant in the P measurements. 

4.2. The effect of cosmic-rays and SBF-measured ct B q estimation 

The HST provides images of exceptional resolution. As a consequence, discriminating 
stars from cosmic rays is a difficult task. This also has implications for SBF measurements. If 
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a cosmic-ray event alters only one pixel, then it is handled as white noise, and consequently, 
it contributes to the Pi signal. But if two or more pixels are affected, then a cosmic ray 
can be confused with a point source and it will therefore contribute to the Pq signal. In the 
WFPC2, the majority of cosmic ray events affect 2-3 pixels, so if undetected cosmic rays 
exist in an image, the SBF results will be affected by an increase of both P\ and Pq. But 
there is noise also due to detected cosmic rays. The pixels containing cosmic rays have been 
masked during the data reduction of the images, and therefore, the noise in these pixels is 
higher than it is in pixels where none of the input images were masked. 

As we will see later, the contribution of cosmic rays to Pq is not negligible, and is 
necessary to control this effect in order to obtain the SBF- measured o"| G . It is for this 
reason that a strategy aimed at estimating and eliminating the contribution of cosmic rays 
from Pq must be developed. Such a strategy is described next. 

In an individual weighted and cosmic-ray cleaned image, Pq has two contributions, faint 
galaxies and cosmic rays: 

P = a 2 BG + al. (9) 

If such individual images averaged, the contributions from galaxies and cosmic rays will 
be different, because of their different natures. Measured in the averaged image, Pq becomes 

^o (av> = + § , (10) 

where N is the number of individual images used for the average. 

In Tables 5 and 6, results for Pq are listed for all the individual weighted and cosmic-ray 
cleaned images and for the averaged images as well. If the influence of cosmic rays were 
negligible, Pq would be the same in all the individual images and in the averaged one. It 
can be seen that this is not the case, from which we may conclude that the contribution of 
cosmic rays to P is non-negligible. 

Equations 9 and 10 can now be used to obtain cr| G and a 2 a . For Pq, the mean of the 
values of the individual images for each chip and filter have been used. Results are listed in 
Table 7 for the two considered values of m c . The final SBF-measured cr| G values for each 
filter can now be obtained from the averages of the single-chip results. They are given in the 
last row of each m c set in Table 7 and will be compared in §5 with the n(m)-estimated <t^ g 
values. 
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4.3. Two consistency tests 

In this section, a couple of consistency tests have been performed in order to check the 
reliability of our SBF-measured cr| G results. 

Comparison of expected and measured Pi 

Together with the PSF-convolved variance, P , SBFs provide the value of Pi, the non-PSF- 
convolved variance. This can be compared with its expected value, directly obtained from 
the read-out noise, the dark current, and the sky brightness analysis of each image. 

For HDF Version 2, each of the weighted and cosmic-ray cleaned images is the result 
of combining several exposures with the same dither position. The different exposures are 
combined with weights proportional to the inverse variance (1/Pi) at the mean background 
level. The variance Pi, in electrons, is computed from the following noise model (Williams 
et al. 1996): 

p 1 = bt + dt + r 2 , (11) 

where t is the exposure time, b is the sky background rate, d is the dark current, and r is 
the read-out noise. The inverse variances, 1/Pi, of each exposure are provided in the header 
of the resulting weighted and cosmic-ray cleaned image. From this information, the value of 
Pi corresponding to the latter can be computed. As an example, these values are listed in 
Table 8 (column 2) for the WF2 images and B450. The Pi values obtained directly from the 
SBF analysis of the images, using m c = 28.8, are listed in column 3. Both computed and 
observed values of Pi are equivalent in all cases. Only a slight excess in the observed Pi is 
noticeable. This excess is produced by cosmic rays, which also contribute to the measured 
P values, as we have shown. 

This test shows that, with this technique, the white noise (Pi) is determined with high 
precision, thereby reinforcing the correctness of Pq measurements. 

Comparison of estimated and measured [o"| G ] 

This second consistency test is based on the comparison of the parameter [o"| G ] computed 
in two different ways. We call [o"| G ] to the variance produced by galaxies with magnitudes 
within a given interval: 
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[a 2 BG ](m A ,m B ) = f* n(f)f 2 df, (12) 
J f B 

where m A is the magnitude corresponding to a flux f A , and to f B - Defined in this 
way, [o"b G ] is the difference between the variances computed with two values of m c , namely 
m c = m A and m c = m B : 

[°bg](™ A,m B ) = al G (m c = m A ) - a BG (m c = m B ) 

(13) 

The n(m)-estimated and SBF-measured [er| G ] can be now compared. To do this, the 
magnitude interval [27.8, 28.8] has been considered. 

As the Williams et al. (1996) photometric catalogue has been used to create the window 
functions, only objects found by them in the magnitude interval [27.8, 28.8] contribute to the 
SBF-measured [cr| G ]. In particular, some of the objects within the interval [27.8, 28.8] in the 
Metcalfe et al. (2001) catalogue will remain masked, namely those which, following Metcalfe 
et al. (2001), are merged in the Williams et al. (1996) catalogue. This implies that the SBF- 
measured [<j| g ] is not expected to coincide with the n(m)-estimated [cr| G ] computed using 
the Metcalfe et al. (2001) data. On the other hand, the n(m)-estimated [crj| G ] computed 
using the Williams et al. (1996) data should, for rigor, be similar to the SBF-measured 
[°bg] onr y if an the objects of the Williams et al. (1996) catalogue in the interval [27.8, 28.8] 
are unmerged. However, if merged galaxies exist in other intervals (Metcalfe et al. 2001), 
there is no reason why they should not be present here also. 

The effect of mergers on [cr BG ] can be tested considering that a fraction of Williams et 
al. (1996) objects in the interval [27.8, 28.8] are the result of a merger between two fainter 
galaxies with integrated fluxes fi and f 2 - We have considered three simple situations: i) 
fi — f2] h) fi = 2/2; and iii) fi = 3f 2 . It should be noted that since SBFs are a measure of 
the second moment of the brightness function, the more similar f 1 and f 2 are, the larger is the 
effect introduced in [cr BG \. So case i) is the most pessimistic and, although it is unrealistic, 
will give the maximum expected effect on [o"| G ] for a given fraction of mergers. 

In Figure 4 we show the results of n(m)-estimated [crg G ] for the three cases and the 
F606W filter, considering different values of the percentage of merged objects. It can be 
seen that, even in the very pessimistic case where 50% of the Williams et al. (1996) objects 
are mergers of two identical galaxies, their influence on [a| G ] is less than 25%. For example, 
in a perhaps more realistic situation in which about 20-25% of the objects are mergers of two 
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galaxies with different magnitudes, the effect on [<r| G ] would smaller than ~10% (see figure 
4). In conclusion, the effect of mergers on [cr| G ] computed using the Williams et al. (1996) 
number counts is small and can be expected to remain below about 15% for any reasonable 
scenario. 

We can hence proceed with our test on [er| G ] for the [27.8, 28.8] interval using the 
Williams et al. (1996) data. The results for SBF-measured and n(m)-estimated [c| G ] are 
given in Table 9. The value of n(m)-estimated [cr| G ] has been computed using Williams et 
al. (1996) and Metcalfe et al. (2001) data. The results are given in columns 2 and 3, while 
column 4 lists the n(m)-measured [o"bg] values. It can be seen that they are not compatible 
with the values obtained from the Metcalfe et al. (2001) data, as expected. 

On the other hand, the n(m)-estimated [<r| G ] obtained from the Williams et al. (1996) 
number counts and the SBF-measured [c^ G ] are very close. In order to compare them with 
more detail, let analyze figure 4 again. The n(m)-measured [cr| G ] value for the F606W filter 
and its error interval have also been plotted in the figure (shadowed region). It can be seen 
now that they fully coincide if a number of mergers about 20-30% is assumed. This situation 
is realistic and compatible with Metcalfe et al. (2001) 's claims. 

Summarizing, the test has been successful and shows that the SBF measurements are 
well-calibrated. 



5. Discussion 

In this section, SBF-measured <r| G results (listed in Table 7) and the n(m)-estimated 
°bg values obtained from both the Williams et al. (1996) and Metcalfe et al. (2001) data 
(listed in Table 3) will be compared. In the following, two possibilities will be considered 
and their consequences discussed: i) that the Williams et al. (1996) data represent the right 
differential number counts, and ii) that the Metcalfe et al. (2001) number counts are correct. 
The obtained SBF measurements will be used to test the validity of these possibilities and, 
as a result, a final n{m) will be proposed. 

5.1. Option A: Assuming Williams et al. (1996) galaxy number counts 

Here we assume that the Williams et al. (1996) number counts are correct. Comparing 
the n(m)-estimated <r| G obtained using Williams et al. (1996) data and SBF-measured 
<t| g , it can be seen that the former are much larger. There are only two possible sources 
to account for this discrepancy: first, a faint unresolved stellar population, belonging to the 
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Milky Way halo, could be responsible of the excess in the SBF signal; and second, the faint 
end of n(m) is different from the fitted one used here to evaluate the n(m)-estimated <t| g . 
In the first case, SBF results may be used to characterize such a halo population. In the 
second case, SBF results may be used to deduce a new faint end of n(m) able to account for 
the SBF-measured <r| G . 

We shall analyze both possibilities in detail and discuss the feasibility of each one and 
its compatibility with the observations; finally, we shall deduce its implications for n(m). 



5.1.1. Faint Milky Way halo stars 

Here, we will consider a Milky Way halo population of faint stars be responsible of the 
observed excess in the SBF signal. We will firstly deduce the halo population necessary to 
cause this SBF signal excess. In order to check the feasibility of this hypothesis, the obtained 
halo population will then be compared with observations in the HDF. 

Lets consider a simple population of objects with absolute magnitude M following the 
standard spatial distribution used by Binney & Tremaine (1987): 

p(r)dr = - ^r^dr, (14) 

w 1 + (r/a) a K ' 

where p(r)dr is the number of objects per pc 3 at a distance from the Milky Way center 
between r and r + dr; a is the core radius, and p is the object density in the Milky Way 
center. For simplicity we take a = 2. 

To derive the SBF signal from the former population, we must first express the equations 
in terms of distance from the Sun (5). This can be done using: 

r 2 = 5 2 + rl - 25r cos6cos/, (15) 

where r is the galactocentric radius, i.e., the distance from the Sun to the Milky Way center, 
and (b, I) are galactic coordinates. The spatial distribution of objects expressed in spherical 
coordinates is then: 

n(5,6,tp)dV 



p a 



2 



a 2 + r 2 + 5 2 — 2r 5 sin 9 cos ip 



S 2 sm9dSd9d(p. (16) 



Integrating for the HDF-N and considering its coordinates / = 127° and b = 54°, the 
former expression is reduced to: 

/ r\ tt 1.38 x 10- 13 p a 2 5 2 
n{6)d6 " a 2 + r 2 + 5 2 - 0.7075r 5 d6 - 
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Now, n(5)d5 can be written in terms of magnitudes by means of the distance modulus 
to obtain: 



n(m)dm 

6.36 x 10- n poa 2 10 a6 ( m - M ) 
" a 2 + r 2 + 100 x \WA™-m) + 7.075r 10°- 2 ( m - M ) 



dm. (18) 



The number of resolved objects with m < 28.8 that should appear in the HDF-N, 
(-^hdf— n^) can ^ e now deduced from eq. 18, as well as the variance that the faint part of 
the population (m > 28.8) would produce (chdf^n^)' 

A7 -(m<28.8) -i 01 w -i n-4 2 

HDF— N = 1-31 x 10 p a 

/28.8-M -i r>0.6x 
dx (19) 
a 2 + r 2 + 100 x 10 - 4 * + 7.075r 10 - 2 * ' 1 ; 



4K 8) = 6.36 x 10- n poa 2 10°- 8(m ^ M) 

roo ^q-0.2x 

X / 28 . 8 „ M a 2 + r 2 + 100 x 10°- 4 - + 7.075r 10°- 2 ^ X ' (20) 

In order to compare these predictions with the HDF data, the V 60 6 filter results have 
been used. A value must be assumed for the core radius, a. Realistic values are around 2000 
pc (Bahcall & Soneira 1980), but 500 pc and 8000 pc have been also used to check a wide 
interval of possibilities. Making ^hdf^n 8 '* ec L u& l to the a^ G excess observed in the filter V 60e 
for m c = 28.8 (that is, 4.87 x 10~ 9 [DN/(s pix)] 2 , see Tables 3 and 7) and introducing the 
value in eq. 20, the central density po can be obtained and used in eqs 14 and 18 to derive 
piocai and N^p 2 ^\ Results are plotted in Fig. 5, where pi OC ai and A^p 28 N 8 ' ) are shown as 
functions of Mf606W, the absolute magnitude of the halo population objects. 

This figure implies the existence of a large number of halo objects that should be present 
in the HDF images. Note that N^p 28 ^ is greater than ~250 in all cases. This result is not 
compatible with the HDF-N observations, where no obvious stars are present, except for a 
few 20th magnitude ones, (Kawaler 1996; Flynn, Gould, & Bahcall 1996). 

In conclusion, the observed excess in SBF-measured a^ G cannot be produced by objects 
belonging to the Milky Way halo. Otherwise a large number of resolved objects from this 
halo population would show up in the HDF-N images, which is not the case. 
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5.1.2. Faint galaxy number counts 

If the observed a BG excess cannot be produced by Milky Way halo objects, the only 
possibility is that it is caused by faint galaxies. The large excess obtained in the SBF- 
measured o~ BG with respect to the n(m)-estimated o~ BG would imply an increase in the slope 
of n{m) at some magnitude fainter than m c = 28.8. This slope can be computed by fitting 
the n(m)-estimated a BG to our SBF-measured a 2 BG and taking the slope as a free parameter. 
If it is assumed that the slope change occurs at m c = 28.8 for all filters, the resulting slopes 
for the fainter range are 7 = 0.60, 0.44, and 0.54 for -B450, V606, and Igu, respectively. These 
slopes would be valid up to fi 450 = 34.4, V 60 6 = 31.9, and J 814 = 32.5 at least since the 
contribution of fainter magnitudes to a BG becomes smaller than the uncertainties in the 
SBF-measured a BG results. If the slope change were to occur at a magnitude fainter than 
28.8, it would result in a steeper n{m). In any case, such big changes in the slope of n(m) 
seem unrealistic. In our opinion, this possibility should be rejected. As a consequence, it 
must be concluded that the Williams et al. (1996) data are incomplete. 

5.2. Option B: Assuming Metcalfe et al. (2001) galaxy number counts 

Assuming that the Metcalfe et al. (2001) differential number counts are correct, the 
SBF-measured cr| G results listed in Table 7 and the n(m)-estimated <r| G values obtained 
using the Metcalfe et al. (2001) data, listed in Table 3, can be compared. 

It can be seen that the SBF-measured and n(m)-estimated <r| G coincide within the error 
bars for the F814W filter, and is very similar for the F450W filter. Only in the filter F606W 
some differences arise. This implies that extrapolation of the Metcalfe et al. (2001) n(m) 
function to magnitudes fainter than 28.8 accounts almost entirely for the measured SBF 
signal, thus indicating a high level of precision in the Metcalfe et al. (2001) data. 

However, as SBF-measured and n(m)-estimated cx| G present slight differences in the 
F450W and F606W filters, the most likely n(m) function that completely fits our SBF 
measurements can be determined. This will be done in the next subsection. 



5.3. The galaxy differential number counts beyond Vq 06 = 28.8 

In this section, the most likely n(m) function will be obtained for magnitudes fainter 
than 28.8. We consider that the most likely n(m) function is that which fits our SBF 
measurements. 
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An n(m) function is completely determined by giving the slope (7) and the number 
of galaxies in a given area and magnitude interval (n m ). In particular we consider n 2 g, 
computed for the magnitude interval [28.75, 29.25] and the area of the HDF-N. Note that 
there are infinite sets of (7, n 2g ) that can produce the same SBF signal. In Figure 6, the 
(7, ^29) pairs that can account for our measured SBF signal have been plotted (solid lines) 
for the F450W, F606W, and F814W filters. With short-dashed lines we represent the (7, 
n 2 g) pairs that would produce the SBF measurements ±ler. In this figure, the (7, n 2 g) sets 
corresponding to the extrapolation of n(m) obtained from Williams et al. (1996) (solid 
circle) and Metcalfe et al. (2001) data (open circle) have also been plotted. 

It can be seen in Figure 6 that the Metcalfe et al. (2001) (7, n 2 g) reproduce our SBF 
results for the F814W filter, as previously shown. For F450W, the Metcalfe et al. (2001) (7, 
n 2 g) is very close to our SBF measurements. In all cases (except for F606W), the Williams 
et al. (1996) (7, n 2 g) results are far from our SBF results, as previously argued. 

The most likely n(m) function can be obtained from Figure 6. We assume that, for each 
filter, the best n{m) estimate for magnitudes fainter than 28.8 is given by the nearest point 
of the solid lines to the Metcalfe et al. (2001) point. In this case, results for the slopes are 
7 = 0.27, 0.21, and 0.26 for B 450 , V 60 6, and J 814 , respectively. 

The results are listed in Table 10 and plotted in Figure 7 (solid lines). The slopes 
obtained are valid down to magnitude 31 at least. The contribution to the SBF signal by 
objects of fainter magnitudes is less than the uncertainty in the SBF measurements. This 
value represents an extension of more than two magnitudes beyond the limits of the previous 
photometric studies by Williams et al. (1996) and Metcalfe et al. (2001). It should also be 
mentioned that, within the uncertainties and except for the objects that could be merged 
into brighter ones in Metcalfe et al. (2001), it is free from incompleteness. 

6. Conclusions 

In this paper, the faint end of the differential galaxy number counts, n(m), has been 
studied by means of SBF measurements. Once the contribution from cosmic rays has been 
evaluated and eliminated from the SBF signal, the background PSF-convolved variance orig- 
inating from faint objects has been carefully analyzed. Our conclusions can be summarized 
as follows: 

• Comparing the SBF-measured <r| G with the n(m)-estimated a^ G predicted by the 
extrapolation of Williams et al. (1996) number counts a clear excess has been found 
in the measured signal. The possibility that the excess might be produced by Milky 
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Way halo stars is ruled out because it would be totally incompatible with the resolved 
stellar population present in the HDF. On the other hand, if this excess is caused by 
a faint galaxy population modifying the faint end of n(m), then the required slopes 
for magnitudes fainter than 28.8 are 7 = 0.60, 0.44, and 0.54 for -B450, Vqqq, and 7 814 , 
respectively. Such big changes in the n(m) slope seem unrealistic. In our opinion, this 
possibility should be rejected. In conclusion, the Williams et al. (1996) number counts 
are not compatible with our SBF measurements, probably owing to the incompleteness 
in their data. 

• Comparing the SBF-measured a^ G with the n(m)-estimated cr| G , predicted by the 
extrapolation of Metcalfe et al. (2001) number counts, we find that they coincide 
within the error bars for the F814W and F450W filters and are similar for F606W. 
This implies that the extrapolation of the Metcalfe et al. (2001) n(m) function to 
magnitudes fainter than 28.8 nearly accounts for the measured SBF signal, indicating 
a high level of precision in the Metcalfe et al. (2001) results. 

• The most likely n(m) function has been obtained fitting our SBF results. Results for 
the n{m) slope for magnitudes fainter than 28.8 are 7 = 0.27, 0.21, and 0.26 for -B450, 
V 60 6, and J 8 i4; respectively. The obtained slopes are valid down to magnitude 31 at 
least. This value represents an extension of more than two magnitudes beyond the 
limits of the previous photometric studies by Williams et al. (1996) and Metcalfe et 
al. (2001). 



This work is based on observations with the NASA/ESA Hubble Space Telescope, ob- 
tained in the Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc. (AURA), under NASA contract NAS5-26555. 
This research has been supported by the Instituto de Astroffsica de Canarias (grant P3/94), 
the DGESIC of the Kingdom of Spain (grant PI97-1438-C02-01), and the DGUI of the 
autonomous government of the Canary Islands (grant PI1999/008). 
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22 24 26 28 30 22 24 26 28 30 22 24 26 28 30 

B 450 AB Magnitude V 606 AB Magnitude I 814 AB Magnitude 

Fig. 1. — Differential number counts of the detected sources by Williams et al. (1996) (filled 
circles) and Metcalfe et al. (2001) (open circles) in the F450W, F606W and F814W filters. 
Solid lines represents the fitted n(m) functions. It can be clearly seen the slope change 
around magnitude 26 found by Williams et al. (1996). 
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Fig. 2. — Example of a SBF measurement: results for the average image in the F450W filter 
with m c = 28.8. Solid line represents the power spectrum fitting. The observed discrepancy 
at low wave numbers between the obtained power spectrum and the fit is due to large 
scale fluctuations in the background brightness of the images. The power-spectrum fitting 
procedure is the following: eq. 2 is used to fit the power spectrum for wave numbers in the 
range [k , /c max ], k nmx being the highest wave number of the computed power spectrum, and 
k a number which we vary from to /c max - As a result, two functions, -Po(^o) and Pi(/co), are 
obtained. The function Po(k ) is also shown (small boxes). It can be seen that this function 
exhibits a "plateau" region. The final adopted result and its uncertainty for P are obtained 
computing the average and standard deviation of -Po(^o) m the plateau interval. For the P\ 
measurement, the procedure is exactly the same as for P . 
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Fig. 3. — P results versus the sky level in WF2 F606W images with m c = 27.8. The absence 
of a trend indicates that flat-fielding errors are not significant in the measurement of Pq. 
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20 40 60 80 100 
Percentage of mergers (%) 

Fig. 4. — Effect on [cbg\ oi a percentage of mergers in the Williams et al. (1996) photometric 
catalogue in the F606W filter. Three simple cases have been considered, i) solid line: fi = f 2 ; 
ii) dashed line: f\ = 2f 2 ; and iii) short-dashed line: fi = 3/ 2 , fi and f 2 being the fluxes of 
the mergered galaxies. The gray region represents the obtained ±1<t SBF-measured [o^g]- 
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Fig. 5. — Local density (pi OC ai) and number of objects in the HDF-N images up to magni- 
tude Vqqq = 28.8 (]Vhdf-n(28.8)) predicted by the considered halo model for three possible 
values of core radius, a. These results have been obtained assuming that the excess in P 
is completely produced by halo objects. It can be seen that the model predictions are not 
compatible with observations. This implies that the observed Po excess cannot be caused by 
halo objects. See text for details. 



-25 - 




2 4 6 8 10 12 2 4 6 8 10 12 2 4 6 8 10 12 

n 29 *10 4 (F450W) n 29 *10 4 (F606W) n 29 *10 4 (F814W) 

Fig. 6. — The n(m) slope, 7, versus the number of galaxies in the HDF-N in the magnitude 
interval [28.75, 29.25], n 29 , for the F450W, F606W, and F814W filters. Solid lines: the pairs 
(7, n 29 ) that would produce our measured SBF signal; short-dashed lines: the pairs (7, n 29 ) 
that would produce the SBF measurements ±lcr; filled circle: the set (7, n 29 ) corresponding 
to the extrapolation of n(m) obtained from the Williams et al. (1996) data; open circle: the 
same corresponding to Metcalfe et al. (2001). It can be seen that the Metcalfe et al. (2001) 
number counts nearly fit out SBF measurements. See text for details. 
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Fig. 7. — n(m) final results in the F450W, F606W and F814W filters. Filled circles represent 
Williams et al. (1996) differential number counts, open circles are the results from Metcalfe 
et al. (2001). Solid line represents the obtained faint end (fainter than magnitude 28.8) of 
n(m). The obtained n(m) function is valid down to magnitude 31, at least. 
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Table 1. The data. 



Image 


Filter 




F450W.dl.dark 


F450W 


17200 


F450W.d2.dark 


F450W 


7900 


F450W.d3.dark 


F450W 


11500 


F450W.d4.dark 


F450W 


12400 


F450W.d5.dark 


F450W 


11800 


F450W.d6.dark 


F450W 


12900 


F450W.d8.dark 


F450W 


13100 


F450W.d9.dark 


F450W 


10700 


F606W.dl.dark 


F606W 


6700 


F606W.d2.dark 


F606W 


4800 


F606W.d3.dark 


F606W 


6450 


F606W.d4.dark 


F606W 


10300 


F606W.d5.dark 


F606W 


15600 


F606W.d6.dark 


F606W 


17300 


F606W.d7.dark 


F606W 


14600 


F606W.d8.dark 


F606W 


10100 


F606W.d9.dark 


F606W 


7100 


F606W.dlO.dark 


F606W 


8300 


F606W.dll.dark 


F606W 


7800 


F814W.dl.dark 


F814W 


10800 


F814W.d2.dark 


F814W 


12200 


F814W.d3.dark 


F814W 


14200 


F814W.d4.dark 


F814W 


12000 


F814W.d5.dark 


F814W 


13000 


F814W.d6.dark 


F814W 


13100 


F814W.d8.dark 


F814W 


5800 


F814W.d9.dark 


F814W 


12900 
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Table 2. HDF photometric zero points. 



Filter 


Chip 


Magnitude (AB system) 


F450W 


PCI 


21.92 


F450W 


WF2 


21.93 


F450W 


WF3 


21.93 


F450W 


WF4 


21.90 


F606W 


PCI 


23.02 


F606W 


WF2 


23.02 


F606W 


WF3 


23.03 


F606W 


WF4 


23.00 


F814W 


PCI 


22.08 


F814W 


WF2 


22.09 


F814W 


WF3 


22.09 


F814W 


WF4 


22.07 
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Table 3. n(m)-estimated cr| G , in units of [ 



m c F450W F606W F814W 





From Williams ct al. (1996) data 




27.8 


2.82 xlCT 9 2.74 xlCT 8 5.50 


xl0~ 9 


28.8 


6.45 xltT 10 6.43 xlCT 9 1.32 


xl0~ 9 




From Metcalfe et al. (2001) data 




27.8 


7.20 xlO" 9 5.76 xl0~ 8 1.05 


xl0~ 8 


28.8 


2.23 xl0~ 9 1.62 xl0~ 8 3.04 


xl0~ 9 
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Table 4. F450W average image (m c = 27.8) SBF results varying the patch radius, in units 



Radius (pixels) 


Po 


Pi 


15 (adopted) 


4.52±0.06 


9.01±0.02 


20 


4.45±0.05 


9.05±0.02 


25 


4.50±0.06 


9.03±0.02 


30 


4.47±0.06 


9.00±0.02 



Table 5. SBF results for m c = 27.8, in units of 1(T 8 x 



Image I PCI WF2 WF3 WF4 

Po Pi Po Pi Po Pi Po Pi 



F450W.dl.dark 


1.90 ± 0.16 


3.86 ± 0.01 


1.05 ± 0.11 


5.4 ± 0.1 


0.86 ± 0.05 


5.4 ± 0.1 


1.09 ± 0.05 


4.6 ± 0.1 


F450W.d2.dark 


3.04 ± 0.10 


7.21 ± 0.01 


1.89 ± 0.07 


11.4 ± 0.1 


1.70 ± 0.23 


10.8± 0.1 


1.13 ± 0.17 


10.0 ± 0.1 


F450W.d3.dark 


1.48 ± 0.06 


5.25 ± 0.01 


1.06 ± 0.10 


8.0 ± 0.1 


1.01 ± 0.07 


7.8 ± 0.1 


1.33 ± 0.10 


6.9 ± 0.1 


F450W.d4.dark 


1.39 ± 0.09 


4.45 ± 0.01 


1.34 ± 0.06 


6.9 ± 0.1 


1.25 ± 0.09 


6.7 ± 0.1 


1.42 ± 0.13 


6.2 ± 0.1 


F450W.d5.dark 


1.73 ± 0.09 


5.18 ± 0.01 


1.21 ± 0.05 


7.5 ± 0.1 


1.36 ± 0.04 


7.2 ± 0.1 


1.89 ± 0.07 


6.3 ± 0.1 


F450W.d6.dark 


1.47 ± 0.06 


4.31 ± 0.01 


1.04 ± 0.12 


6.9 ± 0.1 


1.22 ± 0.03 


6.7 ± 0.1 


1.00 ± 0.08 


6.1 ± 0.1 


F450W.d8.dark 


1.24 ± 0.05 


4.32 ± 0.01 


0.77 ± 0.08 


6.9 ± 0.1 


1.40 ± 0.05 


6.7 ± 0.1 


0.76 ± 0.11 


6.1 ± 0.1 


F450W.d9.dark 


2.10 ± 0.08 


5.95 ± 0.01 


1.13 ± 0.11 


8.9 ± 0.1 


0.98 ± 0.21 


8.5 ± 0.1 


1.70 ± 0.06 


7.7 ± 0.1 


Average image 


0.29 ± 0.03 


0.636 ± 0.001 


0.452± 0.006 


0.901± 0.002 


0.563 ± 0.014 


0.847± 0.004 


0.441 ± 0.006 


0.762 ± 0.002 


F606W.dl.dark 


4.50 ± 0.24 


17.6 ± 0.1 


7.7 ± 0.3 


36.8 ± 0.1 


6.8 ± 0.9 


37.4 ± 0.2 


2.5 ± 0.2 


35.2 ± 0.1 


F606W.d2.dark 


6.1 ± 0.5 


27.2 ± 0.1 


8.2 ± 0.2 


54.9 ± 0.4 


5.5 ± 0.6 


56.5 ± 0.1 


1.4 ± 0.9 


52.5 ± 0.3 


F606W.d3.dark 


4.7 ± 0.6 


22.2 ± 0.1 


5.5 ± 0.5 


42.7 ± 0.2 


6.2 ± 1.7 


42.1 ± 0.3 


8.8 ± 0.6 


42.6 ± 0.3 


F606W.d4.dark 


2.92 ± 0.21 


13.3 ± 0.1 


4.2 ± 0.4 


26.1 ± 0.1 


6.2 ± 0.5 


25.7 ± 0.1 


5.2 ± 0.4 


23.9 ± 0.1 


F606W.d5.dark 


4.71 ± 0.18 


9.3 ± 0.1 


6.1 ± 0.2 


17.6 ± 0.1 


6.8 ± 0.6 


19.1 ± 0.2 


4.5 ± 0.6 


16.8 ± 0.2 


F606W.d6.dark 


2.48 ± 0.09 


7.5 ± 0.1 


6.5 ± 0.3 


16.2 ± 0.1 


6.0 ± 0.2 


17.4 ± 0.1 


5.2 ± 0.4 


14.8 ± 0.1 


F606W.d7.dark 


3.52 ± 0.16 


9.0 ± 0.1 


6.3 ± 0.3 


18.2 ± 0.1 


7.1 ± 0.4 


18.6 ± 0.1 


4.3 ± 0.2 


16.7 ± 0.1 


F606W.d8.dark 


3.6 ± 0.4 


12.6 ± 0.1 


7.4 ± 0.5 


25.2 ± 0.1 


6.7 ± 0.6 


25.7 ± 0.1 


5.6 ± 0.4 


24.2 ± 0.2 


F606W.d9.dark 


5.17 ± 0.24 


20.3 ±0.1 


6.6 ± 0.4 


37.6 ± 0.1 


7.8 ± 0.3 


38.4 ± 0.2 


6.5 ± 0.4 


33.7 ± 0.2 


F606W.dlO.dark 


4.50 ± 0.13 


15.1 ± 0.1 


7.3 ± 0.4 


30.2 ± 0.1 


7.6 ± 0.6 


31.3 ± 0.2 


4.7 ± 0.7 


28.2 ± 0.2 


F606W.dll.dark 


5.37 ± 0.16 


15.5 ± 0.1 


5.5 ± 0.8 


31.3 ± 0.2 


7.4 ± 0.8 


32.3 ± 0.2 


4.4 ± 0.7 


28.9 ± 0.2 


Average image 


0.96 ± 0.18 


1.38 ± 0.01 


3.19 ± 0.13 


2.68 ± 0.02 


3.08 ± 0.20 


2.52 ± 0.04 


3.2 ± 0.5 


1.95 ± 0.14 


F814W.dl.dark 


2.24 ± 0.10 


7.94 ± 0.01 


0.83 ± 0.16 


14.9 ± 0.1 


1.5 ± 0.3 


15.7 ± 0.1 


1.8 ± 0.3 


14.7 ± 0.1 


F814W.d2.dark 


1.89 ± 0.11 


7.59 ± 0.01 


1.99 ± 0.15 


13.3 ± 0.1 


3.23 ± 0.09 


24.4 ± 0.1 


1.70 ± 0.18 


13.3 ± 0.1 


F814W.d3.dark 


2.04 ± 0.04 


6.13 ± 0.01 


1.34 ± 0.17 


11.8 ± 0.1 


1.68 ± 0.11 


11.9 ± 0.1 


1.56 ± 0.07 


10.5 ± 0.1 


F814W.d4.dark 


2.03 ± 0.05 


6.92 ± 0.01 


2.51 ± 0.23 


13.0 ± 0.1 


1.94 ± 0.20 


13.8 ± 0.1 


0.75 ± 0.10 


13.6 ± 0.1 


F814W.d5.dark 


2.1 ± 0.3 


7.00 ± 0.01 


2.6 ± 0.3 


12.9 ± 0.1 


0.92 ± 0.14 


14.6 ± 0.1 


2.07 ± 0.09 


13.0 ± 0.1 


F814W.d6.dark 


2.04 ± 0.11 


6.62 ± 0.01 


2.7 ± 0.3 


12.1 ± 0.1 


3.09 ± 0.17 


13.6 ± 0.1 


1.29 ± 0.16 


12.8 ± 0.1 


F814W.d8.dark 


5.98 ± 0.14 


20.26 ± 0.01 


4.1 ± 0.3 


27.5 ± 0.1 


3.5 ± 0.3 


28.6 ± 0.1 


3.0 ± 0.6 


25.3 ± 0.1 


F814W.d9.dark 


2.55 ± 0.19 


6.77 ± 0.01 


2.0 ± 0.4 


12.6 ± 0.1 


2.3 ± 0.3 


13.0 ± 0.1 


2.60 ± 0.07 


12.0 ± 0.1 


Average image 


0.413 ± 0.019 


1.074 ± 0.002 


1.011 ± 0.024 


1.669 ± 0.007 


0.47 ± 0.05 


2.057 ± 0.012 


0.88 ± 0.03 


1.592± 0.011 



Table 6. SBF results for m c = 28.8, in units of 1(T 8 x [-^l 2 



Image I PCI WF2 WF3 WF4 

Po Pi Po Pi Po Pi Po Pi 



F450W.dl.dark 


1.53 ± 0.04 


3.66 ± 0.01 


0.75 ± 0.08 


5.4 ± 0.1 


0.62 ± 0.06 


5.4 ± 0.1 


0.83 ± 0.04 


4.7 ± 0.1 


F450W.d2.dark 


3.00 ± 0.09 


7.14 ± 0.01 


2.01 ± 0.06 


11.2 ± 0.1 


1.48 ± 0.23 


10.8 ± 0.1 


1.46 ± 0.17 


9.9 ± 0.1 


F450W.d3.dark 


1.43 ± 0.04 


5.24 ± 0.01 


1.04 ± 0.11 


8.0 ± 0.1 


0.81 ± 0.07 


7.8 ± 0.1 


1.22 ± 0.05 


6.9 ± 0.1 


F450W.d4.dark 


1.30 ± 0.06 


4.45 ± 0.01 


1.09 ± 0.06 


6.9 ± 0.1 


1.14 ± 0.05 


6.7 ± 0.1 


1.23 ± 0.10 


6.1 ± 0.1 


F450W.d5.dark 


1.68 ± 0.08 


5.09 ± 0.01 


0.76 ± 0.06 


7.6 ± 0.1 


1.2 ± 0.3 


7.2 ± 0.1 


1.49 ± 0.04 


6.4 ± 0.1 


F450W.d6.dark 


1.46 ± 0.05 


4.31 ± 0.01 


0.59 ± 0.04 


6.9 ± 0.1 


0.67 ± 0.18 


6.8 ± 0.1 


0.77 ± 0.06 


6.1 ± 0.1 


F450W.d8.dark 


1.55 ± 0.04 


4.29 ± 0.01 


0.71 ± 0.07 


6.9 ± 0.1 


1.16 ± 0.07 


6.8 ± 0.1 


0.81 ± 0.07 


6.0 ± 0.1 


F450W.d9.dark 


1.97 ± 0.06 


5.96 ± 0.01 


1.32 ± 0.03 


8.8 ± 0.1 


1.2 ± 0.4 


8.4 ± 0.1 


1.57 ± 0.05 


7.7 ± 0.1 


Average image 


0.265 ± 0.009 


0.626 ± 0.001 


0.317± 0.009 


0.929 ± 0.003 


0.341± 0.015 


0.896± 0.004 


0.247 ± 0.011 


0.822 ± 0.004 


F606W.dl.dark 


4.75 ± 0.17 


17.7 ± 0.1 


5.8 ± 0.3 


37.2 ± 0.1 


6.0 ± 0.6 


37.7 ± 0.2 


2.4 ± 0.7 


35.0 ± 0.3 


F606W.d2.dark 


5.9 ± 0.8 


27.3 ± 0.1 


8.9 ± 0.8 


55.2 ± 0.3 


4.2 ± 0.7 


56.8 ± 0.2 


3.6 ± 0.7 


51.9 ± 0.3 


F606W.d3.dark 


4.9 ± 0.4 


21.2 ± 0.1 


3.8 ± 0.4 


43.0 ± 0.1 


4.7 ± 1.7 


42.5 ± 0.5 


6.0 ± 0.8 


44.8 ± 0.3 


F606W.d4.dark 


2.49 ± 0.22 


13.4 ± 0.1 


1.9 ± 0.4 


26.5 ± 0.1 


4.4 ± 0.5 


26.1 ± 0.2 


4.2 ± 0.4 


24.2 ± 0.1 


F606W.d5.dark 


4.69 ± 0.17 


9.3 ± 0.1 


4.0 ± 0.2 


18.0 ± 0.1 


4.4 ± 0.3 


19.1 ± 0.1 


3.2 ± 0.4 


16.9 ± 0.2 


F606W.d6.dark 


2.41 ± 0.12 


7.6 ± 0.1 


4.2 ± 0.2 


16.6 ± 0.1 


4.1 ± 0.2 


17.6 ± 0.4 


3.8 ± 0.4 


15.1 ± 0.2 


F606W.d7.dark 


2.9 ± 0.4 


9.1 ± 0.1 


4.4 ± 0.3 


18.6 ± 0.1 


4.5 ± 0.3 


19.2 ± 0.1 


3.8 ± 0.7 


16.8 ± 0.3 


F606W.d8.dark 


3.48 ± 0.16 


12.6 ± 0.1 


4.9 ± 0.5 


25.8 ± 0.2 


5.5 ± 1.0 


26.1 ± 0.3 


2.5 ± 1.3 


24.3 ± 0.4 


F606W.d9.dark 


5.3 ± 0.3 


20.6 ± 0.1 


4.9 ± 0.6 


37.6 ± 0.2 


4.4 ± 0.4 


39.2 ± 0.1 


4.8 ± 0.4 


34.1 ± 0.2 


F606W.dlO.dark 


4.41 ± 0.20 


15.0 ± 0.1 


6.2 ± 0.3 


30.5 ± 0.1 


6.4 ± 0.8 


31.7 ± 0.2 


2.6 ± 0.7 


28.7 ± 0.3 


F606W.dll. dark 


5.6 ± 0.4 


15.4 ± 0.1 


3.4 ± 1.0 


31.8 ± 0.3 


5.3 ± 0.4 


33.1 ± 0.1 


2.3 ± 0.5 


29.2 ± 0.2 


Average image 


0.557± 0.021 


1.41 ± 0.01 


1.38 ± 0.15 


2.81 ± 0.02 


1.72 ± 0.10 


2.71 ± 0.02 


1.4 ± 0.3 


2.43 ± 0.09 


F814W.dl.dark 


2.34 ± 0.08 


7.99 ± 0.01 


0.28 ± 0.12 


15.1 ± 0.1 


1.3 ± 0.3 


15.9 ± 0.1 


2.02 ± 0.23 


14.8 ± 0.1 


F814W.d2.dark 


1.51 ± 0.09 


7.64 ± 0.01 


1.26 ± 0.13 


13.3 ± 0.1 


3.6 ± 0.5 


26.2 ± 0.1 


2.09 ± 0.17 


13.2 ± 0.1 


F814W.d3.dark 


2.04 ± 0.04 


6.10 ± 0.01 


1.43 ± 0.09 


11.9 ± 0.1 


1.54 ± 0.22 


12.0 ± 0.1 


1.00 ± 0.09 


10.6 ± 0.1 


F814W.d4.dark 


2.14 ± 0.05 


6.91 ± 0.01 


1.48 ± 0.08 


13.3 ± 0.1 


1.54 ± 0.11 


13.9 ± 0.1 


0.9 ± 0.4 


13.9 ± 0.1 


F814W.d5.dark 


2.1 ± 0.3 


7.02 ± 0.01 


1.4 ± 0.5 


13.1 ± 0.1 


1.37 ± 0.16 


14.9 ± 0.1 


1.6 ± 0.3 


13.2 ± 0.1 


F814W.d6.dark 


1.79 ± 0.06 


6.67 ± 0.01 


2.12 ± 0.12 


12.5 ± 0.1 


2.83 ± 0.15 


13.2 ± 0.1 


1.95 ± 0.12 


12.9 ± 0.1 


F814W.d8.dark 


5.2 ± 0.5 


20.81 ± 0.01 


4.27 ± 0.15 


27.7 ± 0.1 


3.2 ± 0.3 


28.6 ± 0.1 


2.3 ± 0.3 


25.4 ± 0.1 


F814W.d9.dark 


2.61 ± 0.18 


6.78 ± 0.01 


2.13 ± 0.12 


12.3 ± 0.1 


2.50 ± 0.06 


13.0 ± 0.1 


2.04 ± 0.07 


12.3 ± 0.1 


Average image 


0.397 ± 0.022 


1.084 ± 0.002 


0.54 ± 0.03 


1.787 ± 0.007 


0.484±0.022 


2.140 ±0.006 


0.43 ± 0.03 


1.763± 0.010 
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Table 7. SBF-measured cr| G and results, in units of 10 8 x [ 







F450W 


F606W 


F814W 


Chip 




-BO 




-lo 




-BO 










m 


c = 27.8 








PCl( a ) 




0.35 ± 0.16 


0.215 ± 0.007 


3.1 ± 0.9 


0.336 ± 0.020 


0.48 ± 0.10 


0.314 ± 0.006 


WF2 




0.347 ± 0.008 


0.105 ± 0.005 


2.26 ± 0.13 


0.331 ± 0.018 


0.83 ± 0.03 


0.176 ± 0.012 


WF3 




0.468 ± 0.015 


0.094 ± 0.006 


2.72 ± 0.20 


0.36 ± 0.03 


0.52 ± 0.05 


0.218 ± 0.013 


WF4 




0.320 ± 0.008 


0.121 ± 0.005 


3.0 ± 0.5 


0.16 ± 0.05 


0.74 ± 0.04 


0.139 ± 0.021 


SBF-measured 


-BO 


0.37 ± 0.04 




2.8 ± 0.3 




0.64 ± 0.03 




m c = 28.8 


PC1(°) 




0.26 ± 0.05 


0.211 ± 0.006 


0.91 ± 0.11 


0.370 ± 0.011 


0.42 ± 0.12 


0.297 ± 0.012 


WF2 




0.214 ± 0.010 


0.102 ± 0.004 


1.04 ± 0.15 


0.342 ± 0.021 


0.36 ± 0.03 


0.179 ± 0.011 


WF3 




0.242 ± 0.019 


0.099 ± 0.011 


1.40 ± 0.10 


0.32 ± 0.03 


0.23 ± 0.03 


0.251 ± 0.013 


WF4 




0.115 ± 0.012 


0.132 ± 0.004 


1.2 ± 0.3 


0.22 ± 0.04 


0.24 ± 0.03 


0.187 ± 0.013 


SBF-measured 


-BO 


0.208 ± 0.014 




1.13 ± 0.09 




0.31 ± 0.03 





Results from the planetary camera have been scaled to the wide field pixel size. 
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Table 8. Computed and observed P 1 for F450W (WF2) and m c = 28.8, in units of 



Image 


Computed Pi 


Observed Pi 


F450W.dl.dark 


5.2 


5.48 ± 0.02 


F450W.d2.dark 


11.5 


11.25 ± 0.01 


F450W.d3.dark 


7.8 


8.02 ± 0.03 


F450W.d4.dark 


6.9 


6.96 ± 0.02 


F450W.d5.dark 


7.3 


7.62 ± 0.02 


F450W.d6.dark 


6.7 


6.99 ± 0.01 


F450W.d8.dark 


6.6 


6.91 ± 0.02 


F450W.d9.dark 


8.3 


8.81 ± 0.01 
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Table 9. n(m)-estimated and SBF-measured [cr| G ], in units of 10 8 x [ 



Filter 


n(ra)-estimatcd 


n(m)-estimated 


SBF-measured 




using Williams et al. (1996) 


using Metcalfe et al. (2001) 




F450W 


0.22 


0.50 


0.16 ± 0.04 


F606W 


2.10 


4.14 


1.6 ± 0.3 


F814W 


0.42 


0.75 


0.33 ± 0.04 
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Table 10. n(m) Results 



Filter 


7 


SBF limiting magnitude 


F450W 


0.27 


31.0 


F606W 


0.21 


30.7 


F814W 


0.26 


30.8 



